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Abstract 

Tin disturbing functions due to the geopotential and Luni-solar attractions 
are linear and bilinear forms in spherical harmonics. Making use of recurrence 
relations for the solid spherical harmonics and their derivatives, recurrence 
formulas are obtained for high degree terms as function of lower degree for 
any term of those disturbing functions and their derivative with respect to 
any element. The equations obtained are very effective when a numerical 
integration of the equations of motion is appropriate. In analytical theories, 
they provide a fast way of obtaining high degree terms starting from initial 
very simple functions. 


*0n leave from the University of Sao Paulo. 
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Recurrence relations for spherical harmonics 

The Associate Legendre Functions of the first kind 

p"' (x) = (l-x 2 ) m/2 (x 2 -l) a 

* dx^ 


( 1 ) 


satisfy well known recurrence relations, namely 

(Ml) P™ + ,(x) = (2t +1) [x P™ (x)+m(1-x 2 ) 1/2 pf (xlHPj., (2) 


and 

XP™ (x) = P^_, (x) + (t-m+l)(l-x 2 ) 1/2 pj" 1 (x) (3) 


where P^x) = 0 for m>&, according to definition (1). 
Consider now the solid spherical harmonics 

Xj = P'J (sin<j>) e ima r £ 

Vj-P^sin*) 

where 

r 2 = x 2 + y 2 + z 2 , 

X + iy = r cos^ e 1a , 
and z = r si nd> . 


(4) 

(5) 


( 6 ) 

(7) 


Since 0 < m < i , for any i, recurrence relations are necessary for 

a) A given m, increasing values of i 

b) A given pair (m=A-.i), increasing values of 



These distinct cases are necessary in order to avoid singularities in the 
equations we are going to obtain. Such a distinction is not necessary in 
(2) and (3) provided the definition pT = 0 (m >Sl) is assumed, 

Thus we shall provide recurrence relations in the triangular scheme 


m 



Consider first the diagonal recurrence 
(a,?,) U+l.A+l) 


One has 


y £+1 - 

Vl " r 


pJJ] (sin 4>) e 


i (&+1 )a 


But from (2), with x = sin «j 


Pjli - (21+1 ) cos * P l 



I 


I 
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so that 

Xj*] = r ul (2M-1) cost ?l e iXw e ia 
or, in view of (6) , 

x£] * (2t+l)(x+1y) X* 

Similarly 

vj+] = (2*+D p(x + ly) v£ . 

Now, let and consider the vertical recurrence 

(£,m) *v (£+l,m) . 

Combining (2) and (3) one finds: 


From (4): 


3 m B _2£+l p m ( v _ _A+1 p m ( « 
£+1 £-m+l * £-m+l K £-r x ' 


Y m _ ,,£+1 D m Q ima 

V 1 " r P £+l e 


and, using (11) , 


v m _ „£+l 

Vl “ r 


2£+l n m 

wi+r sin ' 1 ' p t 


£+m p m 
£-m+l r £-l 


,ima 


2s.+l .. ..... v m £-Hii .. u m 

r smt a - — r X, 


5,-m+l 


£ £-m+l £-1 


or, using (7) 


m 2£+1 „m _ £+m „ Y m 

' £+1 £-m+l z A £ " £-m+l r A £-l 


Similarly 


£+m 1 v m 

£-m+l r^ y £-1 


(9) 


(10) 


(ID 


(12) 


Y m - ^£+1 Z y 111 
£+1 " £-m+l r 2 T £ 


( 13 ) 


5 


Equations (9), (10), (12), (13), are the necessary relations to be used later. 


Recurrence relations for the cartesian derivatives of Spherical Harmonics 
The solid spherical harmonics can be defined by (Courant-Hilbert, 1973) 


x m B r 2a+l (-0* 

* U-m)! 


L. + i L 
ax ayl 


wr i 




r 

v / 


(14) 


and 


« (-D 

A (t-m)! 


ft- + 1 i-f 

V 

t-m 

r 

,ax a y] 

az, 


r 

, . 


(15) 


From (14): 


£+1 


ytiH-1 _ 2*+3 (-1) 

w " r w. 


l.+llJ 

m+l 

V| l-m 

T 

.ax ay, 


,3zj 

r, 


= r 2ll+3 (-1) 


a 

(-1)* 

’2_ « 1-' 


r j 

T 

ayj 

U-m)! 

.ax ay, 

lazj 


r. 


= r 2t+3 (-1) 


— +i — r " 2 *" 1 x m 

(3x ayj % 


= -r 


*-+i 1 
3X gyj 


Xj + (2t+l) (x + iy) X 1 


*-+1 


ax 


ay) 


= " J 7 x U+i + "7 (2a+l)(x+iy) X nl 


1 v m+l , 1 


Using (14) again: 

,m-l _ 2 k, +3 


ym- I _ t)iTJ ( - 1) 

V 1 ” r 


4+1 


( 4-m+2 ) 


— +i 

a 1 

m-1 

r \ 

a 

t-m ,2 

T 

[ax 

ayj 


,3Zj 

dz^ i 

Tj 


( 16 ) 
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But one has 


o 2 1 

i ? ; 


7 


+ 3 


x „ 2 

T " "3 
r r 


2 2 
3^: 


= _ 

f 2 2 

a . 3 
7T + 7T 

t \ 

1 


i-+l i- 

( 

r • 
1 


ax ay 

r. 


ax ay ( 

,3x ay 

.r 


so that 


y m-l _ 2 ji +3 
*A+1 “ r 


= r 2a+3 


Hi 




fe-H *- 


^mr 3 ^-m m 


= r‘ 


(i-m+2) ! 

(ax 

3yj 

(ax ayj (azj (rj 

( a-m) ! 

r a 

-lS-1 

r -2A-l x m s 

U-rm-2) ! 

[ax 

3yJ 

s. 

1 

(s. 


) v m (2a+l)(x-iy)Xp 


(A-m+2)(j,-m+1 ) 


3x 9yj 


(a-m+2)(jHiH-l) 


D 

3 X 


■i i- 
ayj 


x m . x m-1 + IZi+U ( x . 1y ) x ' 


(17) 


Adding and subtracting (16) and (17), one has 


^ X s! = “ « + U_rn+2)(A “ m+1) X 1+l + 

+ (2£+l) f 0 (18) 

r * 

^ X * = ^7 X £+1 +^r U-rm-2)U-rrtf-l) xj;] + 


+ Zp (2£+l) ^ 


(19) 
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Equations (18) and (19) give the derivatives with respect to x and y. 
In order to find the z-derivative consider again Eq. (14) with &+1 : 


v m 2z+3 

x *+l • r 


U+l -m) ! 


[5- « *J 

m 

f 

a 

ay, 


.32 


9Z 



_ r 2a,+3 (-l)(&-m) 1 3_ r -2A-l x m _ 

(fc+l-m)! az 1 

B _ r 2 _J 8X5L + zm x m 

i-rn+1 az fc-m+1 1 


v 1 " ~ _ ^-rn+1 vfn + z , x m 

az * X)l+1 ^ ' ? * 

In order to obtain the derivatives for the hamonics, one operates 

the X™ and find 


— v™ » - \ y£] + £ U-m+l)U-m+2) v”;] 
ax 

i y™^ + |(*-,n + l Hum*) Y”;] 



U-m+1) Y' 


m 


£+1 


(20) 
as for 

( 21 ) 

( 22 ) 

(23) 
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Disturbi ng function for geopotential and derivatives 
It is given by 

l pa £ 

R ~ l l -r|r Py (sin <{>) [cj cos m x+ S m sin m x] (24) 

£22 m=o r 1 £ 

where 

X = a - 0 

a = right ascension 
0 = Greenwich right ascension 
ij) = declination 

and the other quantities have the usual meaning. 

The cartesian coordinates satisfy the relations 

x + iy = r coseji e 1tx 

z = r sin <j> (25) 

x - iy = r cos <j> e~ 1a 

Let 

«?»»*♦ 5 .i» -o 

so that 

r = i i<. 

z<2 m=0 

Considering that X =■ a-0, one finds 


9 


where 


C? cos m \ + S 1 ! 1 sin in * = (C 1 ? cos m o- s'J sin m o) cos m a + 


+ (cj sin m 0 + S'^ 1 cos mo) sin m a 
= A™ cos m a + sin m a 


r.m 


A^ = cos m e - S™ sin m 0 
B*J = C 1 ! 1 sin m o + S 1 ? sin m o 

A- /V X> 


Also, let us define , from Eq. (5) 


f « uf + 1 V™ 


Thus 


R ” * "»e (A s u t + B i v ™> 


(25) 


(27) 


(28) 


Terms R? can now be generated in succession starting from R° ■ jj- 


H. 

r ’ 


by the diagonal and vertical relations (10) and (13), 
From (10) one finds 


and 


I &+1 

J U1 ■ 


,&+] _ 


24+1 


(X uj - y Vj) 


(29) 




2£+l 
r 


(x V* + y Uj) 


Where, using orbital coordinates, one has 

x = r cost}) cosa = r [cos(w+f) cosn - cos I sin (w+ f) sinn] 
and (30) 

y = r cos<}. cosa = r [cos(w+f) sinn + cos I sin(w+f) cosn] 
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Combining (29) and (28), one has the diagonal recurrence 


0 £+l _ £+1 2£+l 

R *+l • lia e 

or, more conveniently, 


* ^ 


$ K - + b JI] K + K>] 


R e.+1 _ ,,-a+I 2£+l 
Vl ‘ 1 a e r 


2£+l [/..*£+! i_ llt ,£+l\ i|£ 1 /..r>£+l „a£+ 1\ llM /oi\ 

"X (xA a + i yf W u * + (xB £+i “ y/ W \ < 31 > 

r u -» 


where x,y are defined by (30), when necessary, and Aj£j, by (26) 

Note that, in general 


and 


l(xA ni + y B n .') = cj [cos(w+f) cos(n-mo) - cos I sin(oi+f) sin(n-nr)]+ 

F ft ft ft 


+ s n [cos (oj+f ) sin(n-me) + cos I sin(w+f) cos (n-mo)](32) 


^(xB'J - yA 1 ?) = -C 1 ! 1 [cos(wff) sin(fi-me) + cos I sin(w+f) cos(fi-mo)]+ 
r ft ft ft 


+ S 1 " [cos(ui+f ) cos (n-me) - cos I sin(w+f) sin(n-ma)](33) 

Al 


The vertical recurrence relations are also easily found. From (13) it follows 
that 




,m 


£+1 

a+1 


2£+l z 


ru *I 1 


*' n,+1 7 | v m 


£+m 1 

~T~ 

£-m+l r 


U 1 


V' 


,m 

£-1 

m 

£-1 j 


(34) 


where we should remember that , = 0 if m>£-l , Considering (28), 


one finds 



1 ] 


\+] 


* uaf 1 (A 1 


m tl tn 


ii IM + w m \ 

4+1 u 4+l u z+l Vl } 


(35) 


where U ! J +1 , V 1 ^ are given by (34). Also, when orbital elements are used, 

z = r sin $ = r sin I sin (w+f). (36) 

Where using cartesian coordinates, the cartesian derivatives of give directly 
the equations of motion, that is, 


.ni 


X = - i4 + + x (37) 

r J ox 

where X is a nonconservative force (X, Y, Z), and similarly for y, z. When 
using a set of elements (i=l,2,3) and 0^(i~l,2,3), let 


w = col ( x, y, z, x, y, z) 
y " cot ( cti , Qg , Otg , 0-| , 3g , 3g) 

and J the state transition matrix 


d = 


ow 



Define 






where 0,1 are the 3x3 null and identity matrices respectively. 
Lagrange's equations are 



(38) 


( 39 ) 


J 


* 


I 




i 


I 
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where 

F « Jr + R - £ (x 2 + y 2 + z 2 ) (40) 

and 

($) = col (X, Y, Z, 0, 0) . (41) 

One sees that given the cartesian derivatives, Lagrange's equations are easily 
written for any set of elements, provided the state transition matrix J is 
known, for a Keplerian orbit. Or else, one can consider that 

dR _ aR: ax. + aR ay_ + 

act. “ ax 3 a i ay a^ az 3c^ 


and similar equations for 3^. Whatever the form used, If one knows the 
cartesian derivatives of R, the equations are easily constructed by straight 
matrix multiplication. 

Using Eqs. (21), (22), and (23) one finds 


a_ 

ax 



1 

T 


,m+l 
A+l 
m+1 
&+1 , 


1 + A 


U-m+l)U-m+2) 



(42) 


ay 



+ ( A-m+1 ) ( A-m+2 ) 


J v a+1 

1 u m_1 

[ U A+1 


(43) 



( JA-m+l ) 


,m 

Vl 

,»m 

1 


a__ 

az 


(44) 


I 


4 


1 


l 


I 


1 
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Thus , 


OR 


,m 


a(x»y,z) 


VK 




9v 




3(x„y,z) 


+ B 1 


in 


dV 1 


m 


3(x,y,z)j 


(45) 


can be written by recurrence using equations (42), (43), and (44) together 
with (29) and (34), For the geopotential, the recurrence scheme is therefore 
complete. 


Disturbing function for Luni -solar perturbations and derivatives 

For lunar or solar attractions the disturbing function is given by 


where 


R* = l 


i>2 m=0 


L. K » 


(46) 


* 6m 1 •||^f P™ (sin*) P™ (sin *') cos m (o-a 1 ) (47) 


where 6 = gravitational constant 

m‘ = mass of the disturbing body 
c m = 1 (m=0) or 2 (m^O) 

4i 1 ,a ,= declination and right ascension of the disturbing body 

r‘ = earth-moon or earth-sun distance. 

r,<j>,a= spherical equatorial coordinates of the satellite 

Let, from Eq. (4), 


( 48 ) 


4 


*0 


14 


so that 

P™ (sin<f>) cos m a 
and 

Z ! J = r 5, p" 1 (sin4>) sin m a. 


It is easily seen that 


Rf = On' 


*m 


U-m) ! 
U+m)! 



u ; 


m 


+ Z? V 


m y.ra 
a \ " 


(49) 


where U^ m and V^ m are the U^, V™ functions defined over the coordinates of 
the disturbing body, that is 


u i" ” p * < s1n +’> C0SB, “' 

and (50) 

V;" 1 “ Pj (sin«j>') sin m a' . 


Again, any term R^ m can be generated by lower degree terms by recurrence. 

In the equations of motion for the satellite, only cartesian derivatives of W? 

Jo 

and Z* are used. Recurrence relations for U^ m , V^ m are the same than for 
the uj, vj. that is, given by equations (29) and (34) by "priming" all 
variables . 

The recurrence relations for W^, are easily obtained from Eqs. (9) and 
(12), by considering definition (48). One finds, from (9) 

w£| - (2 <l+1)(x w£ - y Zj) 
z£] - (2t+1)(x 
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where x, y are given by (30) when orbital elements are used. Also, from 

( 12 ), 


fw”J 

W H+1 

_ 2)1+1 z 

f „ i 

uHl 

\+1 

- £+m . r 


ylTI 

4 )t+l 

£-m+l 

z m 

r^ij 

&*mi+l 

7 m 

l Vlj 


where z is given by (36) when orbital elements are used. 

The cartesian derivatives of W^,Z^ are found from (18), (19), and (20), 
that is. 


3X 


s, 

-,m 


2~7 


■Vi 1 

,m+l 

H+l 


2 ? 


( 2,-m+2) ( A-m+1 ) 


u-n 

w a+l 

7 * 11-1 

Vl 


. * 

7 


(2A+1) 


yin 

r*J 


(53) 





f 7 m+l 

3 

_ i 

Vl 

sy 

k 

- — 1 
2r 

W m+ ] 

IVIJ 


Ajr (A-m+2) (Ji-m+l ) 
2r^ 


7 m-l 

~h+] 

w m-l 
\+ 1 


+ -^r (2)1+1) 




(54) 


3_ 

3Z 


,m 


,m 


_ - ft-m+1 


Vl 

,m 

^+1 


^ (24+1) 


V 

-»m 


(55) 
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In the above equations the recurrence relations (51) and (52) should be used 
in order to obtain all the derivatives, up to any degree and order, by 
recurrence. 

Therefore, 


aR^ 01 ( (Ji-m)l 

a(x,y,z) = G m> e m [i+m ) : 



3'(x,y,z) 


+ V 


.m 


3Z 111 1 

i 

s{x,y,z) 


(56) 


can be written be recurrence using Equations (53), (54), (55) together with 
(51) and (52) for the satellite terms and using Eqs. (29) and (36) for the lunar 
and solar terms. Obviously, the disturbing function (49) can be obtained by 
recurrence using Eqs. (51), (52), (29), (34). 


Calculation procedure for geopotential 
a) Disturbing function 


Find 


Find 



r \. 

R 3> 



r ° 3 , 

R 2> 

R 3 5 

R 4 > 

R 3' 

R 2 

2 

R 5’ 


by (35) 
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b) Derivatives 


D Rq 

a(x,y,z) 




Find gradient of R 1 * 1 using (45) plus (42), (43), (44), plus (29), (34). 

Xt 


Calculation procedure for luni-solar 
a) Disturbing function 

r'° = ML 


Find R*] , R ' g , R'^, .... by (51) and (29). 
Find 


,3 


I 0 D • Q D fO p I 0 

2 » ’' 3 * n 4 » •••• 


‘ ] » >' 
R*]» R' 


2‘ 

,2 

3’ 


1 R' 1 

2 * ^ ^ ■ 


R' 2 , R' 2 


R 1 

) i • i • 


by (52) and (34). 


b) Derivatives 


8 R' 


TOT ■ [0 '°’ 0] 


Find gradient of R' 1 * 1 using (56) plus (53), (54), (55), plus (51), (52), (29), 

Xf 

(34), for recurrence. 


General comments 

If numerical integration is to be used, the best equations are those 
in cartesian coordinates. The right-hand numbers of the equations 
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+ 


JR , 3R g , 3R° 
ax ax ax 


etc,, can be computed numerically and sequentially by recurrence, Note that 
repeated application of the formulas developed for the cartesian derivatives 
allow partial derivatives up to any order to be quickly computed by recurrence. 
These might be necessary either for construction of variational equations or 
for the integrating process to be used if it requires values of derivatives 
up to a certain order, as in a Runge-Kutta method. 

If an analytical theory is to be developed, the recurrence relations 
can be used with great advantage to generate both the disturbing function 
and its cartesian derivatives by recurrence. If a computer algebraic formula 
processor is available then the right hand member to very high order in analytic 
form can be evaluated very quickly and stored in tape for further manipulation. 
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